Cosmological evolution of interacting phantom (quintessence) 
model in Loop Quantum Gravity 

Puxun Wu and Shuang Nan Zhang ^'^'"^ 
^Department of Physics and Tsinghua Center for Astrophysics, 
Tsinghua University, Beijing 100084, China 
"^School of Sciences and Institute of Math-Physics, 
Central South University of Forestry and Technology, Changsha, Hunan 4IOOO4, China 

^ Key Laboratory of Particle Astrophysics, 
Institute of High Energy Physics, Chinese Academy of Sciences, 
P.O. Box 918-3, Beijmg 1000^9, China 
^ Physics Department, University of Alabama in Huntsville, Huntsville, AL 35899, USA 

Abstract 

The dynamics of interacting dark energy model in loop quantum cosmology (LQC) is studied in 
this paper. The dark energy has a constant equation of state Wx and interacts with dark matter 
through a form 3cH{px + Pm)- We find for quintessence model [wx > — 1) the cosmological 
evolution in LQC is the same as that in classical Einstein cosmology; whereas for phantom dark 
energy {wx < —1), although there are the same critical points in LQC and classical Einstein 
cosmology, loop quantum effect reduces significantly the parameter spacetime {c,Wx) required 
by stability. If parameters c and Wx satisfy the conditions that the critical points are existent 
and stable, the universe will enter an era dominated by dark energy and dark matter with a 
constant energy ratio between them, and accelerate forever; otherwise it will enter an oscillatory 
regime. Comparing our results with the observations we find at la confidence level the universe 
will accelerate forever. 
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I. INTRODUCTION 



Many cosmological observations show that our universe is undergoing an accelerating 
expansion and now mainly consists of two dark components: dark matter and dark energy. 
The dark matter is a clumpy component that traces the baryonic matter and accounts for 
about 23% of present total cosmic energy; the dark energy is an exotic energy with negative 



pressure and accounts for about 72% of tota 



. cosmic energy today. The simplest candidate 



of dark energy is the cosmological constant ij], however it suffers from two problems. One 
is the cosmological constant problem: why is the inferred value of cosmological constant 
so tiny (120 orders of magnitude lower) compared to the typical vacuum energy values 
predicted by the quantum field theory? The other is the coincidence problem: why 
is the dark energy density comparable to the matter density right now? Therefore a 
dynamical scalar field: quintessence 21] is proposed as an alternative of dark energy, but 
it can not explain the region of the equation of state less than —1, which is favored 
by observations 3]. Later Caldwell proposed a phantom field to explain the present 
cosmic accelerating expansion. This field possesses of a negative kinetic energy and so 
has a super negative equation of state. In the Einstein gravity it is found that if the 
universe is dominated by the phantom energy, it will end with a big rip, i.e., a future 



singularity j^. Many works have been done trying to avoid this singularity [6|]. There are 
many other scalar field models: such as quintom [Tl] and hessence [81] . However these scale 
field dark energy models still suffer from the coincidence problem. A possible alleviation 
for this problem is to assume the existence of an interaction between dark matter and 
dark energy j^. 

Recent investigations have shown that there are some new nice features appearing 
in Loop Quantum Cosmology (LQC) (see 
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inflation 



111 ] for recent reviews), such as: easier 



121 ] and correspondence between LQC and braneworld cosmology [13]. The 



LQC is the application in the cosmology context of the Loop Quantum Gravity (LQC) 
(see [l^ for recent reviews), which is a theory trying to quantize the gravity with a non- 
perturbative and background-independent method. Due to the loop quantum effect the 
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standard Friedmann equation can be modified by adding a correction term 
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where H is the Hubble parameter, SttG = 1, p is the total cosmic energy density, pc = 
iQTr^AG'^n, denotes the critical loop quantum density and 7 is the dimensionless Barbero- 
Immirzi parameter (it is suggested that 7 = 0.2375 by the black hole thermodynamics in 
LQG [l8|]). Since this modified equation is correct under the condition that the quantum 
state is semiclassical, this condition is assumed to be satisfied forever in this paper. In 
addition we assume the quantum correlations do not build up during long-term evolution 
of cosmology, otherwise there are additional correction terms from LQC which become 
important [19]. The correction term appearing in Eq. (1) essentially encodes the discrete 
quantum geometric nature of spacetime. When this correction term becomes dominant, 
the universe begins to bounce and then expands backwards. By studying the early universe 
infiation and the fate of future singularity in LQC, it is found that the big bang singularity. 
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16l . 1201 ] . Recently Samart 



the big rip and other future singularities can be avoided [ll|, 
and Gumjudpai 20], Wei and Zhang studied the dynamics of phantom, quintom 
and hessence dark energy models in LQC, and found the results are different from that 
obtained in classical Einstein cosmology. In this paper we will investigate the evolution 
of our universe dominated by a scalar field in LQC, which has constant equation of state 
and interacts with dark matter, and then investigate whether there are some interesting 
features arising from the loop quantum gravity effect. 



II. THE INTERACTING MODEL 



We consider a spatially flat universe in which there are only dark matter and dark en- 
ergy with a constant equation of state w^. Apparently > —1 corresponds a quintessence 
model and < —1 is a phantom case. In addition we assume that between the dark 
matter and dark energy there is an interaction term F. Thus the conservation equations 
for dark matter and dark energy can be expressed as 

+ 3H{1 + w^)p^ = -V , (2) 



3 



Pm + S-ffpm = r , 



(3) 



where p^ and pm correspond to the energy densities of dark energy and dark matter re- 
spectively, and a dot denotes the derivative with respect to cosmic time t. The interacting 
term V is assumed to be F = SHc^p^ + pm), where c is a couphng constant denoting the 
transfer strength. A positive c corresponds to energy transferred from dark energy to 
dark matter and the other way around for a negative one. In this paper we constrain our 
discussion in the case of c > 0. This type of interaction, motivated by analogy with dis- 
sipation of cosmological fluids, has been introduced to solve the coincidence problem [Oj, 



22j, phantom [23| and the (gener- 



and has been studied in the context of quintessence 

I 

alized) Chaplygin gas model [24]. In addition the observational constraints for this type 
interaction dark energy model have been studied in Refs. 25|, l26l |. 

In LQC, using the conservation equation of cosmic total energy p + 3H {p + p) = 0, 
where p = Px + Pm, one can easily obtain the effective modified Raychaudhuri equation 



2 ^' ^ ' V P 
where p is the total pressure (p = WxPx in this paper). 
To analyze the dynamical system, we set 



^ = -4(P + P)fl-2-), (4) 



iV = lna, « = ^, . = ^, (5) 



where Oq = 1 is assumed. Using Eqs.([2l [3], Hj), one can obtain 



u 



-3u{l + Wx) — 3c('U + v) + ?)u{u + WxU + t") ( — 1 H ) , (6) 



2 

v' = —3v + 3c(u + v) + 3v(u + WxU + v){ — 1 ^ ) , (7) 

u + V 



where the prime denotes a derivative with respect to N . As discussed in Refs. jo], [23I, l26| 
this interacting model can solve, or at least ameliorate the coincidence problem in classical 
Einstein cosmological since in the dynamical system there is a late time attractor solution 
with a constant energy ratio between dark energy and dark matter. Therefore, regardless 
the initial conditions, the universe evolves to a final state characterized by a constant dark 
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matter to dark energy ratio. Here we will discuss in LQC whether the dynamics system 
of interacting model exists the attractor solution, and then study the cosmic evolutions 
within different conditions. In order to obtain the possible attractor for the system given 
by Eqs. (jEl [?]), we should firstly solve these equations with u' = and w' = to get the 
critical points: 



Point A : u 



1 I I 4c 

2 2V 



Ur, 



(8) 



Point B : Ur 



--- I — 

2 2V ^ w' 



Ur 



(9) 



Both two critical points correspond to the era dominated by dark matter and dark energy 
with a constant energy ratio between them and exist for c < ^j^. Apparently these 
critical points are the same as that obtained in classical Einstein cosmology [26]. If the 
critical point is stable, it is a late time attractor; otherwise the solution is oscillatory. In 
order to investigate the stability of the critical point, we linearize the system near the 
critical point and arrive at 

Uc + 2f c 



5u' 



3c + 3(— 1 + 2uc + Vc) + + QwxUc 1 



{Uc + Vcf 



Su 



3(c + Uc) + 



[Uc + Vc)"^^ 



6v 



(10) 



6v' 



2vc 

3c + 3(1 + Wx)(-Vc H — 

Uc + Vc 

3 + 3c — 3(1 + Wx)uc — 6vc + 
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Vc{Uc + W^Uc + Vc) 



{Uc + Vcf 



6u 



6v 



(11) 



{uc + ^^c)^ 

Apparently there are two eigenvalues of the coefficient matrix of the above equations. If 
the real parts of two eigenvalues for a critical point are all negative, this critical point 
is stable and is an attractor; otherwise it is unstable and thus oscillatory. We find the 
point B is always unstable, however the point A is an attractor if the equation of state 
for dark energy Wx and the coupling factor c satisfy the conditions l^p^ < c < and 
Wx > —2. Comparing our results with that obtained in Einstein cosmology where point 
A is stable only under condition c < ^j^l j26|, we find that for quintessence dark energy 



Wx > —1 the results in LQC are the same as that in Einstein cosmology if a positive c is 
considered since ^^7^ < 0. However for phantom dark energy Wx < —1 the conclusions 
seem to be different: in region c < ^^p^ or < —2, the point A is stable in Einstein 
cosmology, but it is unstable in LQC, that is, the quantum correction effect will break 
the stability of point A if the interaction factor c is smaller than ^^p^ or the equation 
of state for phantom is less than —2. In Fig. ([T]) we show the stability regions of {c,Wx) 
parameter space. Regions I+II are allowed for Einstein cosmology; however in LQC only 
region II is allowed to obtain a stable solution. 

Since in LQC the interacting quintessence model has the same dynamics as that in 
Einstein cosmology, thereafter we will only discuss the case of phantom < —1- In 
Figs. ([21, [3]), we plot the numerical results for c and satisfying the conditions ^^"^ < 
c < -j^ and Wx > —2. Fig. [2] shows the evolutionary properties of the universe controlled 
by the interacting phantom energy with w = —1.2, c = 0.2 and different initial conditions. 
Apparently the trajectories converge to the same final state determined by parameter c 
and Wx- Since in the final state Qx = Uc and Vim = 1 — Uc, our universe will contain 
both dark matter and phantom energy, and the energy ratio between them approaches a 
constant. Fig. (jS]) shows the evolutionary curve of the equation of state for total cosmic 
energy w = Px/{px + Pm) with Wx = —1.2 and c = 0.2, we find in the final state the 
equation of state is a constant and w > —1, which means that the total energy density 
decreases with the cosmic expansion but the universe accelerates forever. Therefore if 
i±B!L < c < and Wx > —2, regardless the initial conditions, the universe will enter 
a final state with a constant energy ratio between dark energy and dark matter and 
accelerate forever. 

In the following we will give the numerical results for the cases c < ^^7^, Wx < —2 and 
c > ^ according to the Eqs. ([2l [3l iD with pc = 1.5. The cases c < and Wx < -2 
are allowed for the stable solution in classical Einstein cosmology but ruled out by loop 
quantum effect. In Figs. (jU [5]) we plot the evolutionary curves of H{t) and p{t) for the 
case c < with Wx = —1.2 and c = 0.1. In Figs (El [7]) we give the results for the 

case Wx < —2 with Wx = —2.5 and c = 0.25. It is easy to find from these figures that 
at beginning the phantom energy density increases with time, which then leads to the 
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increase of total cosmic energy density. When the total energy density equals to Pc/2, H 
takes the maximum value. When p reaches its maximum value pmax ^ pc, H = and 
then the universe undergos contraction until bounce. Therefore the universe will oscillate 
forever. 

Figs (IHl E]) show the results for the case c > with = —1.2 and c = 0.35 which 
corresponds to the case that the critical points do not exist. Comparing these figures with 
Figs, dll [5], O [7]) , we find, although the universe finally also enters an oscillating regime, 
the process is different from that obtained with the c < or Wx < —2. It is found that 
in this case the energy densities of dark energy and dark matter have the same evolution 
with time, and the H changes from positive to negative (or inverse) when p = ot p pc, 
while in case c < ^-^^ or < —2, H changes from positive to negative (or inverse) only 
at p ^ Pc 



III. CONCLUSION 



In this paper the cosmological evolution with the interacting phantom or quintessence 
dark energy in loop quantum cosmology is studied. We consider the case of dark energy 
with a constant equation of state Wx and the interaction term with the form F = 3Hc{px + 
pm)- It is found that in LQC the dynamic of interacting quintessence model is the same as 
that obtained in classical Einstein cosmological; whereas for interacting phantom model, 
the loop quantum effect reduces significantly the parameter space, in which the attractor 
solution exists. In LQC we obtain the critical point is existent and stable under the 
conditions ^-^j^ < c < and Wx > —2; however in classical Einstein cosmology only 
the condition c < is required. If the coupling parameter c satisfies the stable and 
existent conditions for stable tracking solution, our universe will enter an ere dominated 
by both dark energy and dark matter with a constant energy ratio between them, and 
accelerate forever, although the total energy density decreases with cosmic expansion, 
otherwise the universe will enter an oscillatory regime. 



Recently using the WMAP3 [27[ data Olivares et al 



261] obtained that at la confidence 



level c < 0.0023. More recently in Ref. 25|] by combining the Gold Sne la, BAO and 
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CMB data the authors found that at la confidence level —0.99 < w < —0.83 and c = 
0.0057l'lo;oo26) which show that at la confidence level our universe will enter a final stable 
state and can not oscillate. Letting c = 0.0057 we find in LQC if —1.006 < w < —1 the 
universe with an interaction between dark matter and dark energy will accelerate forever; 
whereas if w < —1.006 it will enter an oscillatory regime. Therefore it is clear that at 
2a confidence level the current observations seem to be unable to predict the late time 
evolution of our universe with the interacting dark energy in LQC . 
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FIG. 1: The stable regions of (wx, c) parameter space. In the region II, critical point A is a 
late time attractor in LQC. In Einstein cosmology critical point A is a late time attractor in the 
region I-|-II. Ill represents the region of the solution without physical meaning. 
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FIG. 2: The phase diagram of the interacting phantom dark energy in LQC with = —1.2 
and c = 0.2. 
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FIG. 3: The evolution of w for total cosmic energy with Wx = —1.2 and c = 0.2. 
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FIG. 5: The evolution of cosmic energy density with time under the condition of c < ^^p^- The 
soHd, dashed and dotted curves correspond to px + Pm, Px and pm respectively. Parameters are 
set as tt; = —1.2, c = 0.1 and pc = 1.5. Obviously px triggers the recollapses, while pm triggers 
the bounces. 
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FIG. 6: The evolution of H with time under the condition of Wx < —2. Parameters are set as 
w = —2.5, c = 0.25 and pc = 1.5. 
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FIG. 7: The evolution of cosmic energy density with time under the condition of c > The 
solid, dashed and dotted curves correspond to px + Pm, Px and pm respectively. Parameters are 
set as w = —2.5, c = 0.25 and pc = 1.5. 
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FIG. 9: The evolution of cosmic energy density with time under the condition of c > The 
soHd, dashed and dotted curves correspond to px + Pm, Px and pm respectively. Parameters are 
set as w = —1.2, c = 0.35 and pc = 1.5. 
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